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Abstract 



In this paper we analyze a singular heat equation of the form $ t + A$ _1 = /. The singular 
term gives rise to very fast diffusion effects. The equation is settled in a smooth bounded 
domain Q C K 3 and complemented with a general dynamic boundary condition of the form 
ai3t — /3Ar$ = where Ar is the Laplace-Beltrami operator and a and j3 are nonnegative 

^vq . coefficients (in particular, the homogeneous Neumann problem given by a = /3 = is included). 

For this problem, we first introduce a suitable weak formulation and prove a related existence 
, result. For more regular initial data, we show that there exists at least one weak solution satisfying 

instantaneous regularization effects which are uniform with respect to the time variable. This 
f ^ ' implies uniqueness for strictly positive times. 

Key words: very fast diffusion, Moser iterations, dynamic boundary conditions. 
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In this paper we address a singular heat equation having the expression 

# t + M~ l = f, (1.1) 

where / is an external force. The equation is settled in a smooth bounded domain ft C M 3 , where the 
restriction to the three-dimensional setting is motivated by physical applications, and is complemented 
with dynamic boundary conditions of the form 

cn?t - pA r $ = <9„iT 1 on T := <9Q, (1.2) 

where Ar is the Laplace-Beltrami operator and a, ft > (in particular they may both be 0, so that 
the homogeneous Neumann problem is included). Equation (jl.l[) can be viewed in the framework of 
nonlinear diffusion problems, i.e., of equations of the form 

■d t - divOT^W) = /, (1.3) 
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with 7Ti £ M. When m > 1, equation (jl.3[) is the well-known porous medium equation |18j (heat 
equation when to = 1); for m < 1, equation (|1.3[) lies in the class of fast diffusion equations (ultra fast 
when 77i < 0); actually, in this regime very fast diffusion occurs in the regions where d is small. In the 
fast diffusion regime, one can also distinguish two sub-regimes. In this context, a crucial role is played 
by the so-called first critical fast diffusion exponent m c := (d is the dimension) which acts as a 
threshold between the good parameter range to £ (m c , 1) and its complementary range to < m c . More 
precisely, when m c < m < 1 and the initial condition $q is non-negative and locally integrable, then 
equation (jl.3p (with / = 0) on the whole space has a unique weak solution which is globally defined, 
positive and smooth. Moreover, solutions emanating from initial data lying in L p {W l ) (p > 1), or even 
in the Marcinkiewicz space M p (R d ), immediately become bounded. The scenario when to < m c is 
drastically different for at least two reasons (we refer to [5] for further remarks). First of all, solutions 
may in general be unbounded and may also be non-smooth. As an example, when f2 = R 3 , one can 
consider, for arbitrary T > 0, the function (see [17] ) 

0(t,x) = 2((T-t)+\x\- 2 ) 1/2 , (1.4) 

which solves (jl.ip with / = and initial condition Qo(x) = 2T 1 / 2 

particular, O is an unbounded function and 0(/j, •) ^ Lf (M. 3 ) for p > 3, at least until it vanishes for 
t > T. In this range of to, the boundedness of the solutions is tied to the summability of the initial 
condition and to the value of to. More precisely, solutions are bounded whenever do £ L p (M. d ) with 
p > p c :— d ^ 2 (see [3]). Note that, when d = 3 and to = — 1 (as in our case (|l.l[l ). p c = 3. A second 
notable feature of the fast diffusion regime, which is already evident in example (|1.4j) . is the possible 
occurrence of extinction in finite time. This means that a solution may become identically zero after 
some finite time T which depends on the initial conditions. Consequently, positivity is lost. The 
ultra fast diffusion regime m < presents further difficulties linked to the mere question of existence 
(see |16j and |llj). In particular, for (nonzero) data in L p with p < p c one may face a phenomenon 
called "immediate extinction" meaning that solutions obtained as limits of reasonable approximation 
schemes can be identically zero for any t > 0. Note that the immediate extinction can occur also for 
boundary value problems with zero Dirichlet conditions. Finally, the case to = to* = jE§, d > 3, 
deserves a particular attention, as observed in the papers [1] and [4] dealing with the asymptotics as 
t /*• T (the extinction time). Note that we always have to* < m c and, for d = 3, we have to* = — 1, 
exactly as in our equation (jl.ll) . 

In this paper, we are interested in the analysis of equation with the dynamic boundary 
condition (|1.2j) in the case when the source term / has zero spatial mean. Then, a straightforward 
computation permits to see that conservation of mass occurs, namely d(t, x) dx = J n do(x) dx for 
any t > 0, which manifestly excludes both immediate extinction and extinction in finite time. Our 
interest is twofold. From the one hand, we aim at proving existence of at least one solution under weak 
conditions on the initial data. To explain what we mean for "weak" , and considering for simplicity the 
case / = 0, we note that, thanks to the monotone structure of (| 1 . ![) - (|1 .2[) . the system admits several 
Liapunov functionals. In particular, physical considerations (see below for more details) suggest to 
introduce an energy functional £ (defined in (|2.12[) below) . Mathematically, the finiteness of £ seems 
to be a minimal regularity condition on the initial datum t9o that allows to define a rigorous concept of 
weak solution (cf. Definition 12.11 below) . Indeed (in the simpler Neumann case a = j3 = 0), finiteness 
of the energy corresponds to asking that i}q — log7?o G L 1 ^). Hence, we have some control on the L 1 - 
norm of the solution (permitting to use Z^-arguments in the analysis) and also a positivity condition. 
Weak solutions emanating from initial data with finite energy £ will be called "energy" solutions. 
Actually, our first result (see Theorem 15 . II below) states that, if the initial data have finite energy and 
/ has zero mean value and satisfies suitable summability conditions, then at least one global energy 
solution exists; moreover, the energy £ ($) remains bounded uniformly in time. 

Once existence is established, we study boundedness and positivity properties of energy so- 
lutions. Assuming the (sole) energy regularity of initial data, we can prove (see Theorem 15. lp that 
t?(i) becomes instantaneously bounded from below, namely, for t > we have 7?(/j) > c(t) > 0, where 
the constant c(t) depends only on the initial energy and may vanish as t — >• 0. If, in addition, i?o 
is (in space dimension d = 3) in some L p -space with p > 3(= p c ), then (see Theorem 15. 4[) we prove 
an analogous bound from above (#(£) < C(i) < +oo for t > 0). These bounds also entail further 
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regularity properties of the solutions, which hold uniformly for large values of the time variable. In- 
deed, c(t) does not vanish, and C(t) does not explode, for t /*• oo. The main tool in our regularity 
proof is a suitably modified Moser iteration scheme. Note that the required extra regularity of the 
initial datum is in complete agreement with the above discussion on the case Q = R d and with the 
explicit solution (|1.4[) . However, since the choice of dynamic boundary conditions Q1.2[) precludes the 
occurrence both of immediate and of finite time extinction, it remains an open challenging question to 
understand whether the restriction p > p c is still optimal in the present setting. Mathematically, the 
main difficulties we encountered in the analysis of this problem come from the choice of the boundary 
conditions (|1.2p which does not allow us to perform some otherwise standard a priori estimates. Ac- 
tually, we have a sort of "asymmetry" of diffusion effects between the equation in the interior domain 
(where the Laplacian acts on and that on the boundary (where the Laplace-Beltrami operator 
acts on (the trace of) •&). Moreover, as the precise statements of our results show, we cannot always 
allow a or /? to be (while they can be simultaneously be since in that case we have the simpler 
case of no- flux conditions). 

A natural application of equation comes from the so-called phase change models of 

Penrose-Fife type [13] • In this physical context, the unknown ■d represents the absolute temperature 
of a material liable to a phase transition, while the source / can also take into account the effects of 
the phase variable on the temperature evolution. More precisely, in the Penrose-Fife model, equation 
(jl.ip is coupled with a parabolic equation of Allen-Calm or Cahn-Hilliard type describing the evolution 
of the phase variable X (see, e.g., [5], [5], [ID], 03], [ED- ^ n this framework, the choice of considering 
a zero-mean-valued forcing function / in (jl.ip can be motivated by the need of replicating in our 
situation some inner cancellation effects that appear in the energy estimates for the full model. An 
application of the present results to the Penrose-Fife system with boundary conditions of the type 
(jl.2p (which has never been studied in the literature, at least up to our knowledge) will be given in a 
forthcoming paper. Besides the Penrose-Fife model, equation (or, more generally, equation ()1.3[) 
with m < 0) comes naturally into play in other physical contexts (see [5] and [7]). For example, it 
appears [T3] in the study of the long-range Van der Waals interactions in thin films that diffuse on a 
solid surface. 

The plan of the paper is as follows. In the next Section [2] we present our assumptions 
and state a rigorous definition of weak solution. Section [3J is devoted to proving the main technical 
lemmas related to the regularization properties of solutions: in particular, we show, by means of 
Moser iteration arguments, that, under suitable assumptions, both -d and satisfy instantaneous 
regularization properties. Section!?] is devoted to showing (global) existence for smooth and bounded 
initial data. Finally, in Section [5] we prove our main results. These regard existence, regularization 
properties, and (in some cases) uniqueness of weak solutions. 

2 Notation and hypotheses 

Let f2 be a smooth bounded domain of R 3 (of course, everything could be easily extended to the one 
and two dimensional cases, where, actually, better results are expected to hold). Let also |f2| = 1 so 
that ||u||ip(n) < |MIl«(o) for all 1 < p < q < +oo, v £ L q (Vl). Let H := L 2 (il), endowed with the 
standard scalar product (■, ■) and norm || • ||. Let also V := H 1 ^). We note by || • ||x the norm in 
the generic Banach space X. and by (-,-}x the duality between X' and X. We will also write || • || p 
for || • ||iP(n) an d || ■ || p ,r for || • ||z>(r), for simplicity. Next, we set Hp := L 2 (T) and Vr := H 1 ^) 
and denote by (■, -)r the scalar product in Hr, by || ■ ||r the corresponding norm, and by (•, -)r the 
duality between Vp and Vr- We also denote by Vr the tangential gradient on T and by Ar the 
Laplace-Beltrami operator. We can thus define the spaces 

n-.= HxH r and V := {z e V : z\ r £ V r } . (2.1) 

Unless specified otherwise, in the sequel we shall make the following convention: when we write h S T-L, 
h will be interpreted as a pair of functions belonging, respectively, to H and to i?r, and both denoted 
by the same letter. On the other hand, when we consider v £ V (or even v £ V), the symbol v will 
be intended, depending on the context, either as a function defined on f2, or as a pair formed by a 
function of fl and its trace on T. 
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For any function, or functional z, denned on ft, we can then set 



ma(z) := / 



n Jn 

where the integral is substituted with the duality (z,l) in case, e.g., z G V . Given a > 0, we also 
define the measure dm, given by 

v dm :— / v + a v, (2-3) 
Jn Jr 

where v represents a generic function in x i 1 (r). Here and below, integrals over T are to be 

intended with respect to the standard surface measure. With some abuse of notation, we will also 
write 

TO ( W ) := foTX iri vdm = -i T mj / udm ' ( 2 - 4 ) 

i.e., the "mean value" of v w.r.t. the measure dm. Here |T| represents the surface measure of T. In 
case a = 0, it is intended that m(v) = toq(w). For p 6 [1, oo) and X a Banach space, we introduce 
the space 

T p (0,+w;X) := \g e Lf oc (0, +oo; X) : sup ||<7lU*( M +i;X) < +00} , (2.5) 

te[o,+oo) ' 

which is a Banach space as it is endowed with the graph norm. Assuming / : (0, +00) x ft — > R be a 
suitable source term and letting a, j3 > 0, we can introduce the system 

t?t - An = /, u = in (0, +00) x ft, (2.6) 

arj t - (3A r n = -d n u, on (0, +00) x T, (2.7) 

0|t=o = #o in ft, (2.8) 

ar)\t=o = &Vo on T. (2-9) 

Our basic assumptions on the initial data are the following: 

j?oGi x (ft), logi?o £ £ 1 (ft), (2.10) 
myoeZ 1 ^), «logr /0 ei'fr). (2.11) 

This very natural condition corresponds to asking that the initial data have Unite energy, where the 
energy functional £ is defined as 

e(#,rj):= [ (tf-logtf)+a / (fy-logT?). (2.12) 
Jn Jr 

In the sequel, we will simply write, with some abuse of language, £ (t), in place of £(iJ(t),r](t)). The 
source term is assumed to satisfy for some given e € (0, 1) 

/ £ L 2 (0,+oo;£ 6 / 5 (ft)) nr 2 (0,+oo ;J L 3+e (ft)), m (/)(i) = 0, for a.a. t G (0, +00). (2.13) 

The reason for taking a source term with zero-mean value is that we need to have a control of the 
spatial average both of $ and of u (and for general / it is not completely clear how to achieve the 
latter). The above assumptions suffice to define a suitable concept of weak solution: 

Definition 2.1. A (global) weak solution (or "energy solution" ) to Problem (P) is a triplet (■&, rj,u) 
satisfying, for all T > 0, the regularity properties 

0eC°([o,T\;L\n)), a V eC°([0,T}- : L\r)) 7 (2.14) 

u E L 2 (0, T; V), u = -1/tf, a.e. in ft, (2.15) 

piog V eL 2 (0,T;V r ), (2.16) 

r) = -l/u r , a.e. on (0,T) x T, (2.17) 
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and fulfilling, for any test function 

£ £ C 1 ([0,T];C°(n))nC ([0,T];C 2 (n)) (2.18) 
and for all times t £ [0, T], the relation 

0(t)£(t) + « / »7(t)£(*) + / / Vu • V£ - / f ?? A r £ 
Jr Jo Jo Jo Jr 

= /V Ht + af I [ [ K+ I f (2-19) 

Jo Jn jo Jr jo j n Jn Jr 

Remark 2.2. It is worth giving some explanation of relation (|2.17p . There, ur denotes the trace 
of m on T, which exists for almost every value of the time variable thanks to the first (|2.15[) . More 
precisely, we have that «r £ L 2 (0, T; i/ 1 / 2 (r)). Actually, for weak solutions we cannot simply write 
rj = ftr since the trace of ft does not necessarily exist. On the other hand, when considering smoother 
solutions (e.g., in the approximation detailed in Section [4] below) it will happen that ft is more regular 
so that rj can be intended, in fact, as the trace of ft. 



3 Main technical lemmas 

In this section we prove some regularization estimates holding for sufficiently smooth solutions of 
Problem (P). In this procedure, a sufficient regularity will always be assumed in such a way that all 
computations we perform make sense. For this reason, we will indicate simply by ft, rather than by 
77, the boundary temperature (indeed, the "bulk" ft has a trace at this regularity level). As a rule, 
we will use the letters c and n to denote generic positive constants, depending only on the set Q and 
on the parameters a and {3, with k used in estimates from below. The values of c, K are allowed 
to vary on occurrence. Finally, Q will denote a generic computable function, increasingly monotone 
with respect to each of its parameters, taking values in [0,+oo). Neither the constants k, c, nor the 
functions Q, are allowed to have an explicit dependence on time. 

3.1 Regularization of u 

Lemma 3.1 (Moser iterations for u). Let S £ R, T > S, a > 0, f3 > 0, and let, for some e £ (0, 1), 



/ £ L 2 {S, T; £ 3+£ (ft)), ||/||L*(s,T;L3+.(n)) = : R ( 3A ) 

Let (u, ft) be a couple of sufficiently smooth functions solving, in a suitable sense, the system 

ft t -Au = f, ft = -1/u, in CI, (3.2) 

aft t - /3A r ft = -d n u, on V, (3.3) 

over the time interval (S, T) . Moreover, let us assume that 

u £ L 2 (S, T; V), \\u\\ La( s,T;V) =■ M, (3.4) 

u(S)eL 1 (n), au(S)EL 1 (T), \\u(S) || L i (n) + a\\u(S) || L i (r) =: U. (3.5) 

Then, for any r £ (0, min{l, T — S}), we have 

IMU~((s+T,T)xfi) < Q{F,M,U,t~ 1 ). (3.6) 



Proof. We will just consider the case when f) = 0, which is more difficult since we cannot get any 
help from the boundary diffusion term. Setting z := —u > 0, we will also assume, for the sake of 
simplicity, that z > 1 almost everywhere. Indeed, if that does not hold, then it is easy to check that 
in the estimates below we can simply replace z with z = max{z, 1}. That said, we test (|3.2j) by z p+1 , 
where p > 1 will be specified later. This gives 
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In order to recover the full V-noim on the left hand side, we can multiply p.7[) by p and then add to 
both hands sides the quantity 

II^IIN NISI- (3-8) 



Hence, using continuity of the embedding V C L 6 (Q) and of the trace operator from V to L 4 (T), we 
get on the left hand side of (|3.7[) a quantity I such that 

1 := k\\ Vz^ f + > K \\z\\l+ 2 +2) + 44 P 2 l p 2 +2) ,r- (3-9) 

We now choose p = 1 for the first iteration. Then, we can estimate the right hand side of (|3.8[) as 

NlSf P =||z||f ^cllzllillzll!. (3.10) 

On the other hand, still for p = 1, we can write 

\f\z 2 < ll/l|3^ 1/2 |||k 3/2 ||6 < a\\z\\ 3 9 + c a \\f\\ 2 M\x, (3.11) 

n 

whence, taking a small enough and collecting (|3.8p - (|3.1ip . (|3.7I) becomes 

^z\\i + a \\ z \\ lx ) + K \\ z \\l + K \\ z \\l r < C (|| Z ||2 + ll/ll^) n^ii^ (3 . 12) 
Thus, integrating over (5*, T) and using Gronwall's lemma, we readily arrive at 

Iklli^syr^CQ)) + a \\ z \\L°°(s,T;LHr)) + lkllL 3 (s\T ; L 3 (n)) + IM|L 3 (s,T ; L 6 (r)) < Q(F, M, U) , (3.13) 

where we can notice that the function Q has no explicit dependence on T. This relation is the starting 
point for the subsequent iterations. 

To proceed, we set r := |i| to be the conjugate exponent of 3 + e. Rewriting (|3.7I) (multiplied 
by p) for a suitable new choice of p, we can now estimate the right hand side as follows: 

P^I/l^^pll/lls+ell^l^pll/lla+ePH^). (3.14) 

Then, adding again the term ||z||(y^2)/2 ^° both hand sides and integrating the result from some t > S 
(to be chosen later) to T, we arrive at 

H Z llL°°(t,T;LP(r!)) + a W Z W P L°°{i,T;LP(T)) + W Z \\ P LP+ 2 {t,T;L 3 (p+ 2 ) + II Z H LP+ 2 (t,T;L 2 (p+ 2 > (r)) 

< c\\z{l)\\l + ca\\z(t)f pr + c(pF + l)PII^ ( 2 +1)(t:T;Lr(p+1)(0)) , (3.15) 
where we used, in particular, that z > 1. Next, we define, for i > and to be chosen later, 

•^t := H^lli^fTi^iLPifn)) H 2; lli,!'i+ 2 (T i ,T,L 3 Pi+ 6 (n)) H z lli,i , i+ 2 (t,T;L 2 (Pi+ 2 )(r))' (3.16) 
Then, elementary interpolation gives 

IMI_L 2 (!>i+i + 1 > ( Ti ,T;L r(p '+i + 1) (Q)) - ll Z llL°°(T i ,T;LPi(I2))ll 2: llLP ! P + 2 (T l ,T,L 3 P. + 6 (0)) 

< /j||2|U°°(Ti,T;£»"i(n)) + (1 - P)ll z lliP I + 2 (n,T,L 3 r'i+ 6 (f2)): (3-17) 

where the index Pi+i and the interpolation exponent p = p(i) are given by the system 

I Pi + 3( Pi +2) r(p 1+1 +l)' 



Dividing the second equation in (|3.18[) by the first one, we have 

1-P \Pi + 2 2 



whence 



Pi 3{pi + 2)Jl-p r 



(3.19) 



and it is easy to compute 
From (|3.20p we also infer 

Being 



<1S±2 = - - i =: K„ (3.20) 
Pi 1 - p r 3 



p ' +2 pz — e (0, 1), provided that Pl > 1. (3.22) 



1 -P = 1 + ^-y ' ( 3 - 23 ) 



we then obtain from the first (|3.18p 



Pi+i = 77- 1 = - H — K e {pi + 2) - 1 



2 1 - p 2 V pi + 2 

where, obviously, H = H(e) > 1 whenever e > 0. 

Taking the pj-th power of (|3 . 1 7[) and using Young's inequality it is not difficult to infer 

Ji > H z llL 2 <Pi+i+ 1 >( Ti ,T;L r(p >+i+ 1) (n))- (3.25) 

We can now start the iteration argument. Let po = 1 and inductively define, for i > 1, Pi+\ := Hpi = 
H z+1 . Moreover, let (for instance), for i > 1, 

3 00 T 

U := t-^, so that V"ti = -. (3.26) 

i=l 

Let us now rewrite (|3.15|) by taking p = Pi+\. Setting also, for brevity, Ji := J Pi and writing t,-|_i in 
place of t (to be chosen below), we then obtain, thanks also to (|3.25p . 

^f+t 1 + "H^ll^c^ +1 ^ : ^^ 1C r» ^ cN^C^^!)!^^- -+- C «||^(^+ 1 )||^:+^ r + C (^ +1 ^ ^ l)^^ 1 ^ 2 . (3.27) 

Let us now make precise the choice of r,. By induction, given Tj, we take Tj+i € (r,, Tj + tj+i) such 
that 



\\z(n + iXl +1 +a\\z(r i+l )\f; i+i)T < -L (||*(*)||£ +1 +a|Kt)||£ +I>r ) 



df 



< c— 
r 



t-j+i 

j2 /-Ti+ti+l / 



" (iK*)iir (Pi+2) +«ik(*)ii?(p 4+a)l r) d * < c V Jf ' (3 - 28) 

where we used that < 2(p* + 2). Thus, ([3~27| - ((3~28| give 

L (r)) 



Jfir + «II*K T<+1 t.l- +1 (D) < « 2H t~" jf + 1 + c(p i+1 F + l) jf + 1+2 . (3.29) 



In particular, we obtain 



^r<c 



(i 2H t- h + + l) J? ' +1+2 . (3.30) 
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Thus, setting m := (H l + 2)/H\ we have 

Ji+i < B? _<i+1 ' Jf +1 , where B { := c{i 2H r- H + H l+1 F + l). (3.31) 
Hence, letting Li := log Ji and Q := logB.;, we can rewrite p.31[) . for i large enough, as 
L l+1 < H-^d + m+1 L t < H-^Ci + Vi+i(H-%-i + ViLi-i) 

i+l i+l i+l 

<---<L l[r 1k +J2{c k - 1 H- k 11 Tfe), (3.32) 
fc=l k=l j=k+l 

where the last productory is understood to be 1 for k = i + 1. A direct check then permits to obtain 
(see also [H]). I 



If u is bounded at the initial time, we can avoid all complications connected with the choice of the 
sequence t;. Indeed, a straightforward modification of the above proof permits to show the following 

Corollary 3.2. Let S G R, T > S, a > 0, /3 > 0, and let, for some e G (0,1), tgTQ hold. Let, as 
before, (u, d) be a couple of sufficiently smooth functions solving, in a suitable sense, system (|3.2[) - (|3.3[) 
over the time interval (S,T). Moreover, let us assume (|3.4p together with 

u(S)€L°°(n), au(S)eL°°(T), \\u(S)\\ L ~ {n) + a\\u(S)\\ L ~ {r) =: U. (3.33) 

Then, we have that 



\\u\\ L ~ ({StT)xQ) <Q(F,M,U). (3.34) 



3.2 Regularization of i3 



Lemma 3.3 (conservation of L p norm). Let S G R, T > S, a > 0, (3 > 0. Let (u, ■&) be a couple of 
sufficiently smooth functions solving in a suitable sense system (|3.2[) - (|3.3[) over Let also, for 

some p > 1, 

f G L 2 (S, T; L 3 (n)), ||/||£ Wi iS(n)) =: F, (3-35) 
tf6l"(S,T i I 1 (n)), W\ L <*>(s,T;mn)) =■ L, (3.36) 
§(S) G LP (fi), atf(S)eL p (r), ||<?(5)IUp(n)+a||i?(5)||w(r)=:e. (3.37) 



Then, 



U°°(s,T;LP(n)) + a||^IU~(s,T;ip(r)) < Q(^, ©)• (3.38) 



Proof. We give the proof only for p > 2. Actually, the case p G [1,2] works similarly but requires 
some change of notation. Then, let us consider first the case when p < 4. We test (|3.2I) by This 
gives 



dt 



11*11? + «ll*H?.r) + II || 2 + /3^^|| V r ^§ || 2 r < p jf (3.39) 



Then, being p < 4, thanks to (|3.36[) . we can add the inequality 

\\#*^\\1<Q{L), (3.40) 

which permits to get the full V-norm of ■d E ?~ from the gradient term. Then, we estimate the right 
hand side of p.39[) as follows: 

P I IJ-I-^^— 1 < || 2 11^?^^ || e < ^ll^ll^c^^^^ ^ c CT || II^H^- (3.41) 
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Taking a small enough, we arrive at the inequality 

+ a \w P ,v) + K\\ni^ 2) < 4mw\ p P + q(l), (3.42) 

whence the thesis follows by integrating in time and using (|3.35[) . p.37p and Gronwall's lemma. In 
the case p > 4, we cannot proceed directly in this way since (|3.36[) does not entail p.40[) . However, 
it is clear that we can first take p = 4 to get an estimate of $ in L°°(S,T; L 4 (il)) (and of ad in 
L°°(S, T; L 4 (r))) and then iterate the procedure finitely many times. I 



Lemma 3.4 (Moscr iterations for z9). Let S G M, T > S, a > 0, (3 > 0. Moreover, if f3 = let 
also a = 0. Let (u, $) be a couple of sufficiently smooth functions solving in a suitable sense system 



(|3.2|> - (|3.3[> over (S,T). Moreover, let us assume that, for some e € (0, 1), 

/ e L 2 (S, T; L 3+£ (tt)), ||/||L»(s,r i £»+.(n)) =: F (3-43) 

G £°°(5, T; L\n)), W\ L co (StTiL i m =: L, (3.44) 

^)£i 3+£ (fi), ^(S)ei 3+£ (r), IIWIU«+«(n) +a||i?(5)||i»+.(r) =: ©• (3-45) 
Then, for any r G (0, min{l, T — S}), we have 

l|0|U~((,sr+r I T)xn) +a|^lli-((s+r,T)xr) < Q(F,L,e,T~ l ). (3.46) 



Proof. To start with a further Moser iteration procedure, we rewrite (|3.39[) . multiply it by p, and 
then add to both sides the term H^ 2 ^ - \\f. This procedure gives 

^ t (w\ p P + <*w\ir) + 4^-2) + p44 p Hp - 2) ,t ^pj n i/i^ -1 + 4n^- (^) 

Then, in the case a = j3 = 0, the proof works similarly with |15[ Lemma 3.5], to which we refer the 
reader for details. In the cases a, (3 > and a = 0, > 0, the argument can be adapted just with 
small variants. Note in particular that, at the first iteration step, we choose p = pa — 3 + e so that 
the latter term in (|3.47[) can be controlled by using Lemma l3~3l I 

In the case when = and a > 0, we can still prove regularization of but the argument is a bit 
more delicate since we cannot take advantage of the boundary gradient term. So, we have to use the 
trace theorem as in Lemma 13.11 This, however, forces us to assume some more summability of the 
initial value of Namely, we have 

Lemma 3.5 (Moscr iterations for #, case a > 0, = 0). Let S G K, T > S, a > 0, and = 0. Let 
(Uj'ff) be a couple of sufficiently smooth functions solving in a suitable sense system l|3.2l) - (|3.3p over 
(S,T). Moreover, let us assume that, for some e G (0, 1), (j3.43p ~ p.44p hold, together with 

^(S)-I 4+£ (J]), tf(S)e£ 4+e (r), \\#(S)\\ L * + c m + WS)\\ L * +c(T) =:e, (3.48) 

Then, for any r G (0, min{l, T - S}), p~46j) holds. 

Proof. We take for simplicity a = 1. The analogue of p.47p reads now 

^(M? + \\n p p ,r) + 4^ lly ^ Pj n l/l^" 1 + c ll*H£=f ■ ( 3 - 49 ) 
Then, we take first p = 4 + e. Owing to Lemma l3~3l we then get 

||0|U-(s,2y^.(n)) + W\l-(s,t;zm*F)) < Q(F, L, 6). (3.50) 

Let us now repeat p.49p for a new p to be chosen below. Let also r = (3 + e')/(2 + e'), where e' is a 
number, also to be chosen below, in the range (0,e]. We then obtain 

pJ \S\^ ^fll/lla+e'll^llr^fll/lla+ell^ll^ 1 .-). (3.51) 
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Hence, integrating p.49j) over (t, T), for some t > S (to be chosen later) and using once more continuity 
of the trace from V to L 4 (T), we arrive (compare with (|3.15[0 at 

ll^lli°°(t,T;£P(n)) + W^W P L^(i,T;LP(r)) + H^llLP- 2 (t,T;L 2 P- 4 (r)) + H^llLP- 2 (t,T;L3p- 6 (0)) 

< ||*(t)||* + Wt)\\ P P ,r + cpF\\n*£_ li(ttTiI *-» m + C ^V- W ^ m 

< \\m\ p P + \\mr P ,r + <pf + mn p L ^ Ht , T;L r< P - 1Hn)r (3.52) 

Now, we perform Moser iterations as in the proof of Lemma 13.11 starting from po = 4 + e. Then, in 
place of (|3.18[) . we get the system 



' 1-p = i 

Pi -2 2(p i+ i-l) ' 

JL _l hl£ = 1 

. Pi ^ 3(p,-2) r(p i+ i-l)- 



(3.53) 



whence one computes, similarly as before, 



and, hnally, 



^-5T^ + 1 -i( 1 + ^5 J,r ')<"-» + i 

= ^^Pi = = ffPi = H ,+ Vo- (3-55) 

Note that, as before, is larger than 1 since e' > 0. 

Then, we can follow with minor variations the proof of Lemma 13.11 up to formula (|3.27[) . In 
particular, Ji is now defined by setting 

Jf' : ~ W^\\ P L^(t,,T:Lp, (O)) + ll^llL°°(T i ,T;Z,Pi(r)) + H^llLP- 2 (n,T;L 2 P ! - 4 (r))~'~ H^llLP- 2 (T i ,T;L 3 Pi- 6 (0))- (3-56) 

Then, the analogue of (|3.28p reads 



ll^(n + i)||^ +1 + H^iJU^.r < 7^ / ' + ll^(t)IIS +1 ,r) dt 

(iw)ii3k- 2 ) + ii*(*)iir (w - 2 ) I r) di < 4 Jf% (3 - 57) 



< c— 

T 



Ti + U + 1 



which holds provided that Pi+\ < 2{pi — 2). As Pi+i is given by (|3.55[) . this means that we need 

(4 + e)iT +1 <2((4 + e)iT-2). (3.58) 
which is easily shown to be true for every i > provided that we choose e' G (0, e] so small that 

(4 + e)H = (4 + e )|±g < 2 ((4 + e) - 2). (3.59) 

Actually, also the choice e' = e works. Notice that, since we used pi < Pi+i < 2(pi — 2), it follows 
that the regularity (|3.48p is needed in order for the iteration scheme to work. From this point on, the 
proof proceeds once again similarly with that of Lemma 13. 11 hence we can omit the details. I 



As before, we have a better result in case the "initial" value of $ is uniformly bounded. The proof 
follows the preceding ones up to straightforward modifications. 

Corollary 3.6. Let S G R, T > S, a > 0, j3 > 0. Let (u, #) be a couple of sufficiently smooth functions 
solving in a suitable sense system (|3.2I) - (|3.3|) over (S,T). Let also (|3.43p - (|3.44p hold, together with 

G a#(S) G L°°(C), ||0(S)IU~(n) + a\\#(S)\\ L ~ (r) =: 9. (3.60) 

Then, we have that 

\\0\\L-«s,T)xn)<Q{L,F,Q). (3.61) 
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3.3 Further regularity of time derivatives 

We now prove that, for smooth solutions of Problem (P), boundedness of $ and u imply some Ir- 
regularity of time derivatives at least for strictly positive times. The proof of this result is a little bit 
tricky in the case j3 > since the diffusion operators on fi and on T act on different functions. In 
particular, we cannot treat the case when /? > and a = because in the derivation of the estimates 
the contribution of the boundary term $ t is explicitly needed. 

Lemma 3.7. Let S G K, T > S, a > 0, /3 > and let d3J} hold for some e G (0, 1). If (3 > 0, then 
let a > 0. Let (u, $) be a couple of sufficiently smooth functions solving in a suitable sense system 
(j3~2>(pT3")) over the time interval (S,T). Moreover, let 

u G L 2 (S, T; V), (3d G L 2 (S, T; V r ), |MU=(s,t ; v) + P\\#\\v(S,Tm < U, (3.62) 

f eL 2 (S,T;H), \\f\\ L 2 {s . T . M) =:F, (3.63) 

u,deL°°((S,T)xn), ||u|U~(( S ,Dxn) + l|0|U-((s,T)xn) ^ U > ( 3 - 64 ) 
for some (given) constants F > 0, U > 0. Then, for any r G (0, min{l, T — S}), we have that 

\\Mz*(S+T,T S H)+<x\\Mz*(S+rF;Hr) + l|u|U-(Sr+r,T;V) + P\\4\\L-(S+r,T ; Vr) < Q{F,U,T- 1 ). (3.65) 

Proof. We just consider the case when a, (3 > 0, which is the most difficult one. For simplicity, 
we can then set a = 13 — 1 and write 



V := / u + / v, (3.66) 
n Jn Jr 

whenever v is, say, an element of H. Then, we test both (|3.2[) and (|3.3[) by u t = fit/"® 2 (note that we 
are always assuming the solution to be smooth enough for our purposes). Then, standard integrations 
by parts lead to 



^ + / r ^WJ £ /« + //£. 



it u ~ ' Jn ^ 2 Jr u ~ Jn 
Of course, using (|3.64|) . the latter term can be simply estimated as follows: 



The first term on the right hand side of (|3.67|) is more delicate. Actually, owing to standard interpo- 
lation inequalities, we have 

< Q(^)ll^llr||Vr^||l r < 0(C/)||^||r||t?|| ff2( r)||^||oc,r 
<Q^)(ll^llr+l|-Ar^||^ + ||^|| 2 



< 



Q(U)(\\0t - A r i?||r-) + Q(U) - Q(U)±\\V r m (3.69) 

Then, the last term is moved to the left hand side and will give a positive contribution. We now 
estimate the first term on the right hand side, which requires some work. Actually, we first notice 
that, comparing terms in equation (|3.3[) . using the trace theorem and elliptic regularity results, we 
have, for arbitrarily small 5 G (0, 1/2), 

Q(U){\\*t - Artfjl*) + Q{U)\\d n u\\ 2 < WIMI? < O(C0ll«ll^+* (n) 

< Q(U)(\\u\\ 2 + ||A«||^_ 4+i(n) + \\u\\ 2 m+S{r) 



< Q(U)(l + || Ak|HM| 2 -*> +\\u\\ 2 H ill ) \\u\\^ (T) ), (3.70) 
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for a suitable p € (0, 2) depending on the choice of S. Then, using (|3.2p and Young's inequality, we 
obtain 



Q(U)\\Aunu\\ 2 -P < Q(U)\\Au\\p < Q(U)\\#t - f\\ p 
<vj^ + QAU){l + \\f\\ 2 )- 



(3.71) 



Moreover, comparing once more terms in (|3.3p and noting in particular that 

A r u = div r (u 2 V r ??), (3.72) 

we arrive at 

||AH| r < Q(U)\\A r nr + Q(t/)||V r tf||ir 

< Q(U)\\A T ti\\ + Q(U) + Q(l7)l|Ar0||r||0||ao,r 

<Q([/)||A r tf||r + Q(t/). (3.73) 
Hence, the last term in (|3.70[) can be controlled this way: 



Q(U)\\u\\ 2 H - 2 ^\\u\\^ 2{r) < Q a (U)\\u\\ 2 Hl{r) +a\\A r u\\ 2 r 



< 



(C0II«IIhi(t) + vQ(U)\\A T n 2 r + aQ(U) 



< Q a {U)\\u\\ 2 m(r) + <tQ(U)\\Mt + cQ(U)\\d n u\\i + aQ(U) 

< QAU)\\n 2 HHr) + vQ(U) j ^ + aQ(U)\\d n u\\l + <rQ(U), (3.74) 



where a > is small and Q a also depends on a. In particular, taking a small enough both in (|3.71[) 
and in (f3T74|) . and collecting ([3J8)) - ([3~74)) . (j3~67| gives 

1 d /„ ll2 f |V r i9| 2 , ,„ „o\ 1 f 



(ll V «l| 2 + / L -^+Q(^)HVr*) 



2dn J r d 2 ^ v /n 1 "V 2 7o^ 2 

<Q(^)(ll/l| 2 +l|tf|lW) + 1 )- (3-75) 
Now, thanks to (|3.62p . for any r as in the statement, we can choose t € (S, S + r) such that 

||V U (t)|| 2 +||Vrtf(t)||^ ^QOAr- 1 ). (3.76) 

The thesis then follows by integrating p.75[) over (t, T) . I 

As before, for regular "initial" data, the property holds starting from the initial time: 

Corollary 3.8. Let S £ R, T > S, a > 0, f3 > and let (J53J hold for some e £ (0, 1). If f3 > 0, then 
let a > 0. Let (u, #) be a couple of sufficiently smooth functions solving in a suitable sense system 
(j3~2>([3~3]) over the time interval (S,T). Moreover, let (|3T52]l - (|3TM|l hoid together with 

u(S)eV, ^(5)GVt, ||u(S)||v + i8||i?(S)||vr=:^ (3-77) 
for some Z > 0. Then, 

\\MV(S,T;H) + a\\0 t \\L><.S,T;Hr) + IMU~(S,T;V) + /^ll^ IU~(S,T;V r ) < Q{P,U,Z). (3.78) 

4 Existence for smooth data 

Based on the previous lemmas, we will prove here that, for sufficiently regular initial data and source 
term, Problem (P) has a global solution in a rather good regularity class. 
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4.1 Local existence of smooth solutions 



We start by showing that, for regular data, a local in time smooth solution exists. For the sake of 
simplicity, we give the proof only in the case when both a and j3 are strictly positive. Actually, the 
other cases can be treated with differences that are mainly of technical character. 

Theorem 4.1. Let a > 0, /3 > 0. Let T > and let 

/eC O ([0,T]xfl), (4.1) 
d eH 2 (fl), fi<Mx) for all x eTi, (4.2) 

for suitable constants < $ < i9. Then, there exists To G (0,T] depending on f and $o (and 
in particular on such that Problem (P) admits a solution (#,it) over the time interval [0,Tq] 

satisfying 

G H 1 (0,T o ;H)nC o ([0,T o \ X IT), = < &(t, x) < 2d for all (t,x) G [0,T ] xfi, (4.3) 
ueL 2 (0,T Q ;H 2 (n)), tfGL 2 (0,T ;ff 2 (r)). (4.4) 

Proof. We only sketch it, since it follows from more or less classical arguments for quasilinear 
parabolic problems. The key idea is to regularize the monotone function 7(7*) = —l/r by introducing 
a function 7^ : R — > R with the following properties 

1R EC 2 (R), 7 fleL°°(R), 7H (r) = 7 (r) Vr G [|,2tf], (4.5) 
3 < m fl < M R such that 771 fl < j' R (r) < M R Vr e R. (4.6) 

In other words, 7i?, coincides with 7 over an interval that is strictly larger than the range of the 
initial values of the problem. Outside that interval, we substitute it by a smooth and bi-Lipschitz 
approximation (e.g., we can take the first order Taylor expansion and then mollify). At this point, 
we can solve the problem by means of a fixed point argument. We just give the highlight of this 
procedure and set, for simplicity, a = /3 = 1. 

As a first step, we take a function c\> in the ball $ of radius 1 of the space L 2 (0, To; H £ (T)), 
where e S (0, 1/2) is fixed, but arbitrary. Then, we solve the initial- value problem 

T] t - Apr/ = 0, rj\ t=0 = # |r, on T. (4.7) 

Thanks to classical results on parabolic equations, the solution 77 is unique and lies in a suitable ball 
of radius Qi(||$o||ff 2 (fi)) 01 the space 

H\0, T ; H e (T)) n L 2 (0, T ; H 2+e (T)). (4.8) 

Note that, indeed, the boundary initial datum lies in H 3 ' 2 (T) by (|4.2[) and the trace theorem. We also 
remark that, by interpolation, the space ()4.8[) is continuously embedded into C°' a ([0, T ] x T) for some 
a > depending on the choice of e (recall that fl is smooth and its boundary T is a 2-dimensional 
manifold). Thus, we can take To so small that, additionally, 

3— 

-= < n(t, x) < -1? for all (t, x) G [0, T ] x T. (4.9) 
As a subsequent step, we solve the problem 

t?t-A 7fl (0)=/, 79| t=0 = t?o, in (4.10) 
7flW = 7«(»»), on T. (4.11) 

This is a quasilinear parabolic system with bi-Lipschitz continuous nonlincarity 7^ and Dirichlet 
boundary conditions. Standard methods permit then to verify that it admits a unique solution & such 
that 

0eff 1 (O,T o ;JT) ) 7nW GL 2 (0,T ;tf 2 (f2)). (4.12) 
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For instance, the basic a priori estimate corresponding to the above regularity can be obtained testing 
the equation in (|4.10[) by the function 9* (7^(1?) —TZ(jn(r]))), where 1Z denotes the harmonic extension 
operator, namely 

- AK(v) =0 in n, TZ(v) =v on T, (4.13) 

where v is, say, a function in Hr- Note that the regularity f|4.8[) of the trace and the C 2 -regularity of 
7# are essential for the sake of obtaining (|4. 12[) . as a direct check permits to verify. 

Then, being 7^ bi-Lipschitz, 7a(i?) lies in a bounded ball of radius Q2(||^o||ff 2 (n)) of the space 

H\0, T ; H) n T 2 (0, T ; H 2 (Q)). (4.14) 

Moreover, the equation in (|4.10p is quasilinear and uniformly parabolic and it has uniformly bounded 
initial data (by (|4.2|l ). boundary data (by f|4.0p ^) . and forcing term (by (|4.ip ). Thus, standard barrier 
arguments entail that, up to possibly taking a smaller initial time To (in a way that only depends on 
the known norms of $0 and /, on the regularized function 7^, and on the truncation values 

I < ${t,x) < 2d for all (t,x) G [0, T ] x O. (4.15) 

Consequently, by (|4.14j) . interpolation, and the trace theorem, the trace function dnj^) lies in a 
bounded ball of radius Q3(||^o||ff 2 (0)) 01 the space 

1/(0, T ; H £ {T)) for some p > 2. (4.16) 

In particular, it is possible to take To so small that the map 

T : $ — !► L 2 (0, T ; H E (T)), T : <j> k> d nl {§), (4.17) 

takes values into Moreover, thanks to (|4.14[) and standard embedding and trace theorems, the map 
T is compact. Finally, a number of standard checks permit to see that it is continuous. Hence, the 
Schauder fixed point theorem applies to T. This gives a local solution to the system 

t ~ A 7fl (tf) = /, in fi, (4.18) 
rjt - A F 7 ? = -d nlR ($), on T, (4.19) 

plus the initial conditions and the boundary condition i?|r = rj. However, thanks to (|4. 1 5[) and the 
latter (|4.5p . 7jj('i9) coincides with 7 (i9) = — 1/# =: u everywhere in [0,To] x O. Hence, (1?, it) is a local 
smooth solution to ()3.2|) - (|3.3[) . This concludes the proof. I 

4.2 Energy estimate 

In this section we prove the basic energy estimate satisfied by solutions of Problem (P). We start by 
recalling a generalized version of Poincare's inequality (see for instance [151 Lemma 3.2] for a proof). 

Lemma 4.2. Assume £1 is a bounded open subset ofM. d . Suppose v £ W 1 ' 1 ^) and v > a.e. in fi. 
Then, setting K := J n (\ogv) + , the following estimate holds: 

\\v\\ LH n) < \n\e ClK + ^\\Vv\\ LHnh (4.20) 
the positive constants C± and C2 depending only on Q. 

Lemma 4.3 (Energy estimate). Let a > 0, p 1 > and let ([2TTU|) - ([2TTT|) and ([233]) hold. Let (i?,it) be 
a sufficiently smooth solution to Problem (P). Then, we have that 

£{t) + [ (|| V^l 2 + 0|| V r \ogf}\\ 2 r ) + \\u\\ 2 L2(tit+1 . v) < Q(Eo), V< > 0. (4.21) 
Jo 



14 



PROOF. Test (j2l)f and (j2~7f with 1 + u. This gives 

±£{t) + {\\Vuf + /3||V r lognl) = J f(l + u). (4.22) 
Then, using (|2.13|) and the Poincarc-Wirtinger inequality, we estimate the right hand side as follows: 
[ f(l + u)= [ fu = [ f(u- mn(u)) < \\\Vu\\ 2 + c||/|| 2 6/5(n) . (4.23) 
Integrating (14.221) in time from to t, we obtain 

£(t) + J (l\\Vu\\ 2 + /3||V r lo g m) < E + c||/||i a(0iT . x6/8(n))) Vt > 0. (4.24) 

To obtain the control on the full V-norm of u. we use Lemma 14.21 Thus, taking v = l/# = — it in 
P~2T))) . we infer 

||u|| £ i (n) < |n| e Ci/„iog-* + ^L|| Vti || il(n) . ( 4. 25) 
Hence, squaring, integrating over (t, t + 1), and using (|4.24[) . we obtain 

\\u\\ L 2 {t , t+1 . L i m < Q(E ). (4.26) 

The thesis follows combining (|4.26|) and (|4.24p . I 

4.3 Proof of global existence 

Our next aim is to prove that, for smooth initial data, the solution constructed in Section [4. II has, in 
fact, a global in time character. We can treat all cases with the exception of a = 0, /3 > 0. This case 
will be dealt with separately in Section 15.31 below. 

Theorem 4.4. Let a > 0, /3 > and, if a = 0, then let ft = 0. Let ([2~T3")) hold together with 



i? eV, g <i?o(x) < for all a; G ft, (4.27) 

for some < $ < i?. Then, there exists a global solution to Problem (P) satisfying, for all T > 0, 

i? G F 1 (0,T;i/)n J L oo ((0,T) x 0) n L°°(0, T; V), (4.28) 

tt G tf^O, T; ff) n L°°((0, T) x SI) fl L°°(0, T; V) n L 2 (0, T; H 3/2 {Q)), (4.29) 

otj G H^lOjT; Hr), fin G L°°(0, T; Vr) H £ 2 (0, T; ff 2 (r)). (4.30) 

Moreover, if either a > and /3>Oora = /3 = 0, then we have more precisely 

G L 2 (0, T; if 2 (ft)), it G L 2 (0, T; i? 2 (ft)). (4.31) 



Proof. Let us fix (an arbitrary) T > 0. In order to get a local solution via Theorem 14. II we need 
to construct sequences of regularized initial and source data {f?o.n} and {f n } satisfying, for all n G N, 
(|4.1j) and, respectively, (|4.2p . Moreover, we need do. n — > an d f n — > f hi proper ways. The details 
of the rcgularization procedure arc sketched in Section [5.21 below, to which we refer the reader. 

Then, thanks to Theorem 14.11 for all n > we have a solution to the n-Problem defined on 
some interval (0, To !Tl ), with To, n < T. We now deduce global in time estimates independent of n and, 
for the sake of simplicity, we shall directly work on the time interval (0,T) rather than on (0,To jn ). 
As usual, this can be justified a posteriori by means of standard extension arguments. 

First of all, we can apply Lemma 14.31 which gives the estimates 

II^I|l»(0,T;L1(S))) + a\\r)n\\L°°(0,T;Li(r)) + ||«n|U=»(0,T;V) + P\\ l°g Vn II L^(0,T;V r ) < c > ( 4 - 32 ) 
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for c > independent of n. Indeed, from this point on, we go back to the notation rj n when we indicate 
the boundary value of d n - Next, we can apply Corollaries 13.21 13761 (with S = 0), which give 



IKi||l°=((0,T)xQ) + ||$n|U°°((0,T)xn) < c - 

Clearly, the same uniform boundedness properties hold also for the traces on T. 

Next, by Corollarv l3.8l (here the restriction on a and j3 comes into play), we have 



\\&n\\m(0,T;H) + a||f?n||i/i(0,T;ffr) + IKi IU°° (0,T;V) + P\\Vn\\L°°(0,T;Vr) < c - 

Comparing terms in equation (|3.2[) . we also obtain 

\\&Un\\L 2 (Q,T;H) < C. 

Moreover, in the case when /3 > (and hence a > 0), we observe that 



u„\r 



and 

Vn 



Vn 



< C, 



(4.33) 

(4.34) 
(4.35) 

(4.36) 



L 2 (0,T-V r ) 



the latter bound following from (|4.32j) - (|4.33[) . Hence, by standard regularity results for elliptic prob- 
lems with Dirichlet boundary conditions (see, e.g., [5J Theorem 3.1.5]), we obtain 



Kll 



L 2 (0,T;H 3 / 2 (n)) 



< c. 



Thanks to (|4.32p and <|4.35[> . we can apply the trace theorem ® Theorem 2.7.7], which yields 

\\9nU n \\L 2 (0,T;H r ) < C 

Consequently, using the regularity of the boundary initial datum (|4.27[) we obtain 

\\Vn,t\\L 2 (0,T;Hr) + l r /« \\l 2 (0,T;H 2 (T)) < c - 

Next, we observe that 

A r = -rAri)n - 2- 



Thus, using the boundary analogue of (|4.33p . (|4.39[) . and the Gagliardo-Nirenberg inequality 

IMIw^<c|!<4 2 ||<^+c|MI^, 

which gives that |Vr?7n| 2 G ^ 2 (r) uniformly w.r.t. n, we readily obtain that 

1 

Vn L 2 (0,T:H 2 {T)) 



< c, 



(4.37) 

(4.38) 

(4.39) 
(4.40) 

(4.41) 

(4.42) 
(4.43) 



whence flj/JH) can be improved to 

\Wn\\L 2 (0.T-H 2 (Q,)) < C 

and the same bound holds for {i) n } thanks again to (|4.33p and (|4.41l) . 

Let us now deal with the case when /3 = 0. Then, if it is also a — 0, we are just dealing with 
no- flux conditions. Hence, from (|4.35[) we directly deduce (|4.43[) and, consequently, (|4.31[) . If, instead, 
/3 = and a > 0, we then have that d n u n = —arj n> t, so that, from (|4.34[) . f|4.35[) . and regularity 
results for Neumann elliptic boundary value problems, we deduce once more (|4.37p . However, it does 
not seem possible to arrive at (|4.43p since we do not have sufficient regularity of d n u n . 

In all cases, standard applications of the Aubin-Lions lemma permit to take the limit n /* oo 
and obtain (up to the extraction of a subsequence) existence of a triplet 77, u) solving Problem (P) 
and complying with the regularity properties (|4.28p - (|4.30p and, possibly, (|4.3ip . In particular, it is 
worth noting that, in the case a = j3 = there is no boundary function r\. Otherwise, r\ coincides 
with the trace of i? thanks to the quoted regularity properties and to continuity of the trace operator. 
The proof is concluded. I 
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5 Weak solutions 



5.1 Existence for finite-energy data 

Theorem 5.1. Let assumptions (|2.10[) - (|2.11[) and (|2.13[) hold and let us assume that, if a = 0, then 
,8 = 0. Then, Problem (P) admits at least one energy solution u) satisfying the uniform energy 
estimate ^ 

£{t)+ l (\\\ Vu \\ 2 + P\\Vr\ogi9\\ 2 r ) <E + c||/||i 2(0it . i6/5(n)) , Vt>0. (5.1) 
and the further regularity property 

\\u(t)\\ LO , m <Q(E ,M tl T~ 1 ) fora.e. t>r, re (0,1), (5.2) 
where Eo represents the "initial energy" 

E Q :=£(i) ,m) (5-3) 

and we have set 

M e ■= ll/llL2(0, + oo;L<5/5(fi)) + ll/!lr 2 (0,+oc s L3+«(O))- (5-4) 

Proof. As before, we start with the case when a and (3 are both strictly positive. The other cases 
can be treated with small variants that will be outlined at the end of the proof. In order to apply 
Theorem 14.41 we consider a sequence of initial data {$o,n} sucn that 

i?o, n £ V, < ± n < i?o,n(z) < < +oo, a.e. in fl and on T, (5.5) 
#o,n->-tfo strongly in L 1 (Sl), $o, n \r -> Vo strongly in L 1 ^), (5.6) 

as n / oo. The construction of such a sequence is sketched in Section I5~2l below. Notice that we do 
not need to approximate /. Then, for any n > 0, the assumptions of Theorem l4.4l are fulfilled and we 
have existence of a "smooth" solution (z9„, u„). Moreover, we can still use the energy estimate, which 
gives the analogue of (|4.32[) . As a further consequence, we can apply Lemma 13. li on the generic time 
interval (t, t+ 1), so that, combining (|4.2ip with (|3.6p . we obtain 

||un(*)||z,°°(n) < Q(Eo,-r _1 ) for a.e. t > t > 0, (5.7) 

where it is worth noting once more that Q does not depend on the final time T . In other words, we 
have a uniform bound. 

At this point, we use -^-techniques in order to take the limit n /* oo. To this aim, we work on 
the generic interval (0,T) and rewrite the approximate equation (|2.6D in the equivalent form (which 
is possible since u n has the good regularity properties (|4.29p ~) 

d t u n = ulAu n + u 2 n f. (5.8) 

Then, we test (JSHJ) with <p £ L°°(0, T; Hq(Q)), where we choose s > § so that H$(fl) C C°(TT), 
continuously. Note that, since ip = on V, we do not have to take into account the contribution of 
the boundary. Using (|4.21[) and (|5.7[) . it is a standard matter to get 

\\Un,t\\LHT,T;H-') < C- (5.9) 

Next, for any e > and r € M. we introduce the approximate sign function as sign e (w) := £ ^„| ■ Then, 
we write the approximate equations (|2.6p and (|2.7I) for the indexes m and n, take the differences and 
test them, respectively, by sign e (u m — u n ) and by sign 6 (u mi r — u n ,r)- Indeed, for any m,n eN and 
e > 0, we have 

sign e (u m - un) e L 2 (0, T; V r ) (5.10) 
since sign e is Lipschitz. Then, we arrive at 

((#m,i - ^ n ,t),sign £ (u m - u n )) + a((n m . t - rj n ,t), sign £ (u m -u n )) r 

- (3(A r (n m - rj n ),sign £ (u m - u n )) T < 0. (5-H) 
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Here, we used integration by parts to prove that 

- (A(u m - u n ), sign £ (u TO - u n )) > -(d n (u m - u„), sign £ (w,„ - u n )) r , (5-12) 

which is possible thanks to the good regularity of approximate solutions, to monotonicity of sign e , 
and to (|5.10[) . Thus, integrating (|5. 10[) over (0,i), for a generic t G (0, T], we obtain 



((0m,t - i?n,t),sign £ (u TO - u n )) + ot I ((r) m>t - r/„ jt ), sign e (u m - u n )) r 

Jo 

-p (Ar(r?m - 7? n ),sign £ (u m -u n )) v < 0. (5.13) 



Next, we take the limit e \, and notice that, since u h> -l/wis monotone increasing, it turns out that 
sign (i? m - #„) = sign (-u m - u n ) almost everywhere in (0, T)xVL and sign (?7 m - rj n ) = sign (u rn - u n ) 
almost everywhere in (0,T) x L, where sign denotes the sign- like function with sign o (0) = 0. Then, 



((^m,t - #n.t), sign (tf m - ■&„)) + a / ((r) m ,t - r/„,t), sign (r? m - r/„)) r 

Jo 

- P / (A r (r/ m - »7n),sign (?7 m - r] n )) r < 0. (5.14) 



Applying the Brezis-Strauss lemma (see [5J Lemma 2]), we can integrate by parts the boundary 
Laplacian, which gives a nonnegative contribution. Thus, we arrive at 

\\& m (t) - ^(OHi + a\\ri m {t) - ?7„(t)|| lr < ||l? ,rn - ^o.n^ + a||^o,m - ^o,re|| ir - (5.15) 

Taking the supremum for t G [0, T], we finally obtain 

i9 n ->-tf strongly in C°([0,T];£ 1 (n)), ar] n -> my strongly in C°([0, T]; X 1 ^)), (5.16) 

for suitable limit functions i? and ?y. Moreover, using (|4.21[) . (|5.9[) . and the Aubin-Lions lemma, we 
have that, up to the extraction of a subsequence, 

u n ^u weakly in L 2 {0,T;V) and strongly in i 2 (r, T; iJ 1 ^^)), (5.17) 

for all r G (0, 1), 8 G (0, 1), and for some limit function u. In particular, taking a test function £ 
satisfying (|2.18[) . we can test both (|2.6[) and (|2.7p (at the n-level) by £ and integrate by parts. This 
gives the n-analogue of (|2.19[) . Moreover, the convergence properties proved above permit to take the 
limit n oo to show that (|2.19[) holds also at the limit level. Hence, to conclude the proof it remains 
to identify the functions i9 and 77 in terms of u. To do this, we first observe that, combining the first 
(15.161) with (|5.17[) . it follows that d = —1/u almost everywhere in (0,T) x f2. 

Next, by (]5 . 1 T[) and continuity of the trace operator, we have, for all r G (0,T), 

K)r -> u r strongly in L 2 (t, T; H 1/2 ~ S {T)), (5.18) 

with Mr denoting the trace of u. In particular, thanks to arbitrariness of r, pointwise convergence 
holds on (0,T) x T, up to the extraction of a further subsequence. Being r/ n = — l/(u n )r almost 
everywhere in (0, T) x T and for all n G N, we then deduce that 

T] n -> — - a.e. on (0,T) x T. (5.19) 

Combining this with the second (|5.16|) . we then obtain relation (|2.17|) . This concludes the proof of 
the Theorem in the case a > 0, P > 0. 

Let us now consider the case a > and P = (the case a = p = is essentially already 
known and, in any case, it is simpler to treat). Since Lemma 13.11 docs not depend on whether P is 
zero or not, estimate (|5.7I) (and consequently (|5.9[) 1 still holds when P = 0. Moreover, it is easy to 
check that the Cauchy argument (|5.10[) - (|5.16[) can be reproduced also when P = 0. The proof of the 
Theorem is concluded. I 
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Remark 5.2. Note that, in the case when a > and j3 > 0, 770 needs not be the trace of i?o- Thus, 
there is a boundary layer in the sense specified by (|5.18[) - (|5.19[) . However, (|2.17[) still make sense 
thanks to instantaneous regularization properties. 

Remark 5.3. The reason for not considering the case where a = and (3 > can be already seen in 
the Definition ^. ll of a weak solution. Indeed, if a = and (3 > 0, then the last integral on the first row 
of (|2.19p may make no sense since t? needs not have any summability property on the boundary. On 
the other hand, if a > 0, then that integral is finite thanks to the latter (|2.14p . Actually, in the case 
when a = and (3 > 0, we can still arrive at formal estimates of the form (|5.1[) and (|5.2[) . However, 
these estimate do not suffice to prove existence exactly for the reason outlined above. 



5.2 Approximation of data 

In this part, we sketch the approximation of data needed in the proofs of Theorems 14.41 and 15.11 Since 
the procedures are rather standard, we just give the highlights without entering too much into details. 
Approximation of data for Theorem 14.41 We start with the initial data. Let us given $q S V 
such that d < $o(x) < $ for all x S S7 (as in the statement, see (|4.27p ). and let 770 be its trace on 
r. Since we want to apply Theorem 14.11 (cf., in particular, (|4.2[1 ), for any n s N we need to have an 
approximate datum 7? , n (with trace 770, n) such that $o in £ H 2 (Q) for all n and 

1. tV„ — > ?? m V, 

2. t? < tV™(aO < # for all x E Ti. 

To construct ?9o,n we then consider the elliptic problem 

U,,n - £A0„,n = 0o in n, 

\-^<9„i\n = -±A r r)o,n + Vo,n - Vo on r > 0o,n|r = r]o,n- 

Then, the standard elliptic theory gives that, for any fixed n, $o,n lies in i/ 2 (SI) and 770, n lies in H 2 (T). 

To prove the convergence $o,n — > $0, w c test (|5.20|) |o by (7?o,n — A#o in ) and, correspondingly, 
(I5.20l) |r by (1 + 7i _1 )?7o,n — Vo- Then, by integrations by parts and Young's inequality, it is not difficult 
to obtain 

||i?o,n|ft + 2n|| f to I „-»jb||r- < WMl, 

which clearly implies the desired convergence. On the other hand, to prove uniform boundcdness of 
1^0, n, we use a maximum principle argument. We prove only the upper bound, the lower one being 
completely equivalent. For a generic M > 7?, we test (|5.20p |o by (^o,n — M) + . Integrating by parts 
the Laplacian and using (|5.20p |r. we then obtain 



l?O, n (0O.n - M) + + / VoAVO.n - M) + < / M\n~M)++ / % (V0,n ~ Af) + • 

: Jr Jn Jr 

Adding to both sides the quantity — M J n (i?o,n — M) + — M f r (i]o.n — M) + , we arrive at 

||(7\„-A/)+|| 2 + ||(77 ,„-M)+||2 < [ (# -M){# , n -M)++ /"(77o-M)(77o,„-M)+<0. (5.21) 



n Jr 

This clearly gives the desired inequality 7? 0i „(a;) < ■& for all x G SI. 

Concerning the source term, recalling that / satisfies (|2.13p , a combination of truncation and 

mollification techniques, together with a suitable correction of the spatial mean values permits to 
construct a family of functions {f n }, n £ N, such that 

/„ e C°([0,T] x O), ||/„|U»(( 0l T)xn) < n, (5.22) 

/„ / strongly in L 1 ((0,T) x SI), (5.23) 

(/n(t))n = for aU t€ [0,T], (5.24) 

||/n|U 2 ((o,r ; i 3 + e (n)) < c. (5.25) 
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The details of this construction are left to the reader. Notice that (|5.22j) - (|5.25[) suffice both to apply 
Theorem 14. II at the level n to get a local smooth solution and to perform the estimates of Theorem l4.4l 
uniformly in n in order to let n — > oc. 

Approximation of data for Theorem 15.11 We detail such an approximation just in the case 
a > 0, which is a little bit trickier. So, let (i?o,%) G ^(fl) x L x (r) such that log7? G L^O) and 
log 770 G i 1 (r). We need to construct t?o.n in such a way that properties (|5.5[) - (|5.6|) hold. Then, we 
consider first the function $q. Using that T is smooth, we can first extend it (e.g., by reflection) 
to a neighbourhood f2* of ft. It is then clear that the new function, note it as f?o,*, lies in L 1 (0*); 
moreover, log^o.* G L 1 ^*). Next, we truncate i?o,*j setting 



^On := mm \ max { $o,*j — \ ,n> \^ j a.e. in (5.26) 



1 

n 

Finally, we regularize, setting i?q 2 ^ := p n * tf^nj where {p n } is a suitable sequence of smooth and 
compactly supported mollificrs. Then, straightforward checks (based on the properties of convolutions 

(2*) 

and on Lebesgue's theorem) and a standard diagonal argument permit to verify that t9 n is smooth 
and tends to $ strongly in L 1 (£7). Moreover, based on Jensen's inequality, it is not difficult to verify 
that, for some c > independent of n, 

log-tfg<c(l + ^log-tf ). (5.27) 

Finally, we pass to the boundary component. First of all, thanks to smoothness of T, we can find 
e > such that 770 can be extended (e.g., constantly along directions orthogonal to Y) to a function 
770.* defined on a neighbourhood r,5 := {x G M 3 : d(x,T) < 5} of T. Thanks to Fubini's theorem, it 
is then clear that both 770,* and log 7/0,* lie in L (Tg). Then, we truncate 770,* (as in (|5.26[) ) obtaining 
Van (which can be seen as a function defined on the whole of R 3 ). Next, we mollify 77^, introducing 
77^ := p n * 1]^, for p n supported, say, on the ball _B(0,1/tt,). Finally, we take a cutoff function 
tpn G C°°(R 3 ; [0, 1]) such that ip n is identically one on T 1 / 2n and ip n is supported on T 1 / n . Then, we 
se t "Holi := Von^Pn in such a way that rj^ n belongs to C°°(IR 3 ) and t/q 3 ^ tends to in L 1 (fi), while its 
trace tends to 770 in i 1 (r). Moreover, as above, one can check that log~(77Q 3 ^) is uniformly controlled 
in -/^(r) in the sense of (|5.27|) . Then, the required approximation of $0 is obtained simply taking 

UQ,n ■— ^o,n ' 'l0,n- 

5.3 Solutions with regularizing effects for 

In this last part, we extend the previous results in three directions. First, we prove that, if t? (and, 
possibly, 770) enjoy higher summability properties, then there exist weak solutions whose component 
1? satisfies time-rcgularization properties in the spirit of Lemmas l3.4l and l3.5l Second, we demonstrate 
that, under the same type of conditions on the initial data, existence holds also for a = and /3 > 
(recall that we could not deal with this case for L 1 initial data, cf. Theorem 15. ip . Third, we see that 
uniqueness holds in the class of solutions with regularizing effects. All these results are collected in 
the following 

Theorem 5.4. Let assumptions (|2.10j) - (|2.11j) and (|2.13|) hold and let in addition, for some e € (0, 1), 
/ e L 2 (0, +00; L 3+e (ft)), N s := ||/|U»(o, +aoi is+.(n)). (5-28) 
Moreover, if a > and j3 = 0, iet 

tfo g L 4+e (n), a7?o G L 4+C (r), (5.29) 

whereas in the other cases let 

7? G£ 3+e (ft), ai?o G L 3+C (r). (5.30) 
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Then, Problem (P) admits at least one energy solution ($,it) satisfying (|5.1[) . (|5.2p . together with the 
regularization estimate 

WmW^m^Q^o+e^^T- 1 ) V*>t, re(0,l), (5.31) 

where we have set 

E 0+e := £(,?„) + ||*o|li£. (n) +a\\M%t.<r>, ( 5 - 32 ) 

the exponents 3 + e being all replaced by 4 + e in the case when a > and /? = 0. Moreover, in ah 
cases with the exception of a = and /? > 0, we have 

||$i|U 2 (t,t+l;#) + Ot\\r) t \\L2(t,t+UH r ) + IMI L~(t, +00; V) + /3|M|£°°(t,oo;V r ) 

< Q(E 0+Cl N e ,T~ 1 ) Vt>r, re (0,1). (5.33) 

Finally, whenever there exists an energy solution satisfying (|5.3ip and (|5.33[) . then it is the unique 
solution in that regularity class. 

Proof. For the sake of simplicity, we just prove the theorem by directly working on the "limit" 
solutions without referring to an explicit approximation scheme. That said, we first observe that the 
energy estimate (|4.2ip still holds. Moreover, we can still rely on the conclusion of Lemma \3. II Next, 
thanks to assumption (|5.28p . we can apply Lemma |3 . 31 over the generic time interval (0, T). This gives 

IW*)IU3 +e(n) +a||0(t)||i»+.(r) < Q(iV e ,E 0+e ) (5.34) 

(here and below, 3 + e is replaced by 4 + e in case a > and /3 = 0). Hence, for t > r > 0, we can 
apply Lemma [33 (or Lemma l3~5|) over the generic time interval (t, t+1). This gives (|5.31[) . 

Finally, in all cases with the exception of a = and j3 > 0, we can apply Lemma 13.71 over the 
generic time interval (t, t + 1) where t > r > 0, which gives (|5.33[) . 

Then, to conclude the proof, it just remain to show that, in the case when a = and /3 > 
(that we set equal to 1 for simplicity), a weak solution still exists under the above assumptions. To 
this aim, we consider the system 

§ n ,t - Au„ = /, d n = -l/u n , in O, (5.35) 

-r) n ,t - A r ??„ = -d n u, on T, (5.36) 
n 

complemented with the usual initial conditions. Then, for all n £ N, there exists at least one weak 
solution (d n ,u n ). Moreover, thanks to (|5.33[) and to regularity arguments similar to those performed 
in Section |4~31 (z9„, u n ) is smooth enough in order for the system to make sense in the above "strong" 
form, at least on time intervals of the form (r, T) for all r > 0. 

In addition to that, we still have estimates (|5 . 1 [) and (|5.2p . Moreover, it is worth noting that 
(|5.2p holds independently of n. Indeed, looking back at the proof of Lemma 13.11 it is immediate to 
check that estimate (|3.6p is independent of a. Actually, when one performs the iteration argument 
(see p.29p ) one simply has a functional Ji that depends on a, but the bound for such a functional 
remains unchanged. The same holds when applying Lemma |3. 41 In conclusion, we have the uniform 
bound 

l|tfn(*)|U»(n)4j«n(*)IU»(fi) ^(Eo+e^r- 1 ) V t > T, 7 6(0,1). (5.37) 

Thus, thanks also to Lemma T3. 31 for any r G (0, 1), T > r, we have 

d n weakly star in L°°(0, T; L 3+e (Q)) and weakly star in i°°((r,T) x fl). (5.38) 

Moreover, since f3 > 0, as an additional consequence of estimate (|3.39p (with p = 3 + e) we have 

1+1 3+e 

Ir" 2 lli 2 (0,T;V) + II 7 ?™ 2 llL 2 (0,T;V r ) - C ' (5.39) 

In particular, being 

W„ = -^-tf^W^, (5.40) 
1 + e 
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we have, from (|5.38[) - (|5.39[) . 



l + e 



||W n || , 3+« <c W„ a r , rnTlff J^n 2 I 6+2, <c. (5.41) 

In particular, this fact tells us that, in the present regularity setting, rj n can be directly seen as the 
trace of "d n . More precisely, applying the trace theorem, we have 

||tm|| 1+^3+, <c. (5.42) 

L 2 (0,T;W3+«' 2 (r)) 



On the other hand, testing (|3.2p by a generic function <f> £ Jf^fi) of unit norm and recalling (|5.1[) 
and (|5.28p . we obtain 

tfn.t-^i weakly in L 2 (0,T; H^iQ)). (5.43) 
Thus, using (|5.4ip . (|5.43|) and the Aubin-Lions lemma, we infer 

■d n strongly in L 2 (0, T; W/ 1- ' 7 ' 1 ^ (ft)) (5.44) 

for all <7 6 (0, 1). Hence, by continuity of the trace operator, 

T] n -^t] = d\r strongly in L 2 (0, T; L P {T)) for some p > 1. (5.45) 

Thus, writing the weak formulation (|2.19[) of system (|5.35[) - (|5.36[) . it is immediate to check that one 
can take the limit n — > oo therein. In particular, it can be standardly proved that 

u n ^u = -1/0 strongly in 77(0, T; L p (ft)), say, for all p E [1, 2). (5.46) 

Hence, in place of (|2.17[) . we can now directly write rj = -$|r in the sense of traces and a.e. in (0, T). 
Moreover, thanks to the additional regularity properties coming from (|5.41[) - (|5.42[) . we could even 
relax a bit the requirements (|2.18[) on the test function £ (we omit the details). 

Finally, we come to uniqueness. Let us given a couple of energy solutions ui) and ($2> 1*2) 
both satisfying (|5.31|) and (|5.33[) . Thanks to smoothness for strictly positive times, we can proceed as 
in Section [5.11 testing the difference of the equations by sign e (ui — 112). We then integrate over (r, T) 
for r > and arrive at the analogue of (|5.13[) . namely 

(0?l,t - i?2,t),sign e (ui -u 2 )) + ol J ((tjij - 7]2j),sign E (ui - u 2 )) r 

-13 j (A r (m - 77 2 ),sign £ ( Wl -«a)) r < 0. (5.47) 

Note that we cannot integrate directly over (0, t) since (|5.3ip and (|5.33p do not extend to r = 0; in 
other words, we do not have sufficient regularity to use sign e (Mi — u 2 ) as a test function over (0,T). 
Hence, we first need to take the limit e \ 0, obtaining 

pi{t) - + HM*) ~ *»(*)Hi,r ^ IN r ) ~ Will + - ^{r)\\ 1T . (5.48) 

Letting r \, and noting that energy solutions are continuous with values in L 1 (cf. (|2.14p ). we obtain 
the assert. The proof is concluded. I 

Remark 5.5. In principle, the uniqueness proof works also in case ((5~30|) (or ([5^29]) ) is not satisfied, 
i.e., under the conditions of Theorem 15. II However, in that case the statement may be vacuous since 
we do not know whether there exist energy solutions satisfying (|5.31l) and (|5.33p . Actually, even in 
case (|5.30p (or (I5.29P - ) holds, the proved properties does not exclude that there might exist other 
energy solutions of Problem (P) that do not regularize with respect to time. The same observation 
can also be referred to the case when a = and (3 > since (|5.33p is not known to hold under these 
conditions. 
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